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Abstract 

We consider single band of conduction electrons interacting with displacements of the transitional 
ions. In the classical regime strong enough coupling transforms the harmonic elastic energy for 
an ion to the one of the well with two deep minima, so that the system is described in terms of 
Ising spins. Inter-site interactions via the exchange by electrons order spins at lower temperatures. 
Extension to the quantum regime is discussed. Below the CDW-transition the energy spectrum of 
electrons remains metallic because the structural vector Q and the FS sizes are not related. Large 
values of the CDW gap seen in the tunneling experiments find their natural explanation as due to 
the deep energy minima in the bound two-well electron-ion complex. 
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I. INTRODUCTION 



Origin of charge density waves (CDW) in the transition-metals_^ dichalcogenides 
(TM2HDC) remains enigmatic since their discovery in the early 70's [l]. At the time, 
the most popular theoretical scenario was the so-called "nesting" , the congruency of two or 
more Fermi surfaces (FS) separated by a vector Q in the momentum space. Interactions 
with the momentum transfer Q would lead to the CDW instability (with Q becoming the 
structural vector) and to the opening of energy gaps on FS. 

The alternative explanation [2^ ascribes the structural instability to the presence of saddle 
points in the electronic spectrum near the Fermi energy. The logarithmic singularities in the 
electronic density of states (DOS) at the saddle points favor instabilities with the momentum 
transfer, Q connecting the two saddle points. The mechanism also leads to the energy gaps 
in the electron spectrum. 

Subsequent experiments and band structure calculations for real materials gave no sup- 
port to either of the two concepts. The CDW transition does not affect properties of 
dichalcogenides noticeable. Materials remain metallic below the transition temperature, 
TcDW- Superconductivity in 2iJ-TaSe2 and 2if-NbSe2 takes place on the background of the 
CDW phase. 

The density functional calculations jsl (2if-NbSe2) have shown that the contribution from 
the nested pieces of FS into the charge susceptibility is negligible at the expected value of 
the vector Q, although, indeed, the ground state of the system at T = is the CDW phase. 
In js, ^ it was hypothesized that CDW transitions in 2i7-dichalcogenides are driven by a 
phonon mode that becomes soft due to strong electron-phonon {e-ph) interactions, provided, 
though, that the e-ph matrix element is enhanced at the experimentally known structural 
vector, Q. 

The ...oscopic phy.c= g expla... ...the. dffle.enc. between p.opert.e. of 2H- 

TaSe2 and 2iJ-NbSe2, nor the CDW gap values ~ O.leV [5[ which are surprisingly large 
compared to the transition temperatures, TcnwiTa) = 122K and TcDw{Nb) = 33. 5K. (In 
the infrared data (2iJ-TaSe2) the gap of such order of magnitude is noticeable even close 
to the transition temperature). 

The mean field analysis of the CDW-transitions performed in 0, O] has shown that in 
2i^-TaSe2 the coherence length is surprisingly short: ^0 ~ 3 — 10 A. (The CDW transition 
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in 2if-NbSe2 demonstrates features more typical for a mean field structural transition). 

We address these peculiarities below. We argue that one may discern two stages at the 
CDW- formation, at least tentatively. At first, strong e-ph interactions locally bind electrons 
and ions together in kind of a polaronic effect which was first discussed for VsSi in [9]. The 
potential in which moves an ion becomes anharmonic. At a strong enough e-ph coupling 
the potential possesses the two deep minima, so that the system can be described in terms 
of Ising spins. The subsequent ordering of the local sites at a lower temperature is due to 
inter-site interactions via the exchange by electrons. "Trapping" of electrons by the elastic 
field is responsible for the main energy gains and determines, whereby, the CDW gap by 
the order of magnitude. The available experiments are discussed below in the light of these 
results. 



II. CHOICE OF THE MODEL. 

Before proceeding further, we briefly summarize the information pertaining to the CDW- 
transitions in TM2HGC. The mean-field parameters governing the transition were discussed 
in 0]. The lattice superstructure below Tcdw [l] is formed by the triple-Q modulations of 
the ionic positions along the three symmetry axes of the hexagonal Sif-phase. The Q-vector 
may be incommensurate, but even then Q remains very close to a*/3 (here a* = An/y/Sa). 

It was perceived already in the early 70'sjl^ that the structural changes in 2if-NbSe2 are 
related to displacements of the niobium ions. Correspondingly, in what follows, we choose a 
simplified model and consider only shifts of the transition-metals ions (Nb or Ta) along each 
of these three symmetry lines. The dramatic softening of the longitudinal acoustic mode 



with Q about a* /3 was recently observed in Note, in passing, that the optical modes 



and the longitudinal acoustic branch of the same Si- symmetry are actually linearly coupled 
near point a*/3 of the Brillouin Zone (BZ). We return to this later to show that softening 
of the optical mode is accompanied by softening of the acoustic phonon, and vice versa. 

Transition-metals-dichalcogenides are layered compounds with the weak Van der Waals 
coupling between the layers. We restrict ourselves by the quasi-two-dimensional properties 
(Q2D) of a single layer. 

We choose a single isotropic band of electrons interacting strongly with displacements of 
the transition-metal ions. Implicitly, the band is assumed to be predominantly of the 4f- 
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(Nb) or 5f- (Ta) - character. 

In the model 3,|4| the "bare" frequency of the propagating phonon mode, oo^^k) is renor- 
mahzed by the polarization operator 12|] (assuming the e-ph matrix element g{p,p + k)is 
peaked at k ^ Q). Unlike js], Q], our attention is concentrated on the local environment of 
a single ion. We assume that the optical phonon frequency has a negligible dispersion at 
k = Q. In this respect, as we discuss later, our model shows some features common to the 



Holstein model 
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For a heavy ion one may expect that quantum effects are not crucial for properties of the 
lattice, at least, at not too low a temperature. The temperature of the (incommensurate) 
transition in 2iJ-TaSe2, Tcdw = 122ii' [l| is rather high and, in the first approximation, 
one can neglect the kinetic energy of the Ta ions. 

With the notations uj for a displacement of the single ion at the point Ri the "bare" 
elastic matrix is: U{Ri — Rj)uiUj. At i = j it is merely the potential energy of an oscillator: 

U{Ri - Rj)uiUj MujIu^J2 = kul/2. (1) 

Interactions with electrons change the elastic matrix. We calculate the contribution to the 
elastic energy that is due to interactions of electrons with arbitrary static displacements. 
With the Hamiltonian for e-ph interaction in the form 

He-ph = Y.gu^d^^{Ri)^{R,) (2) 

i 

the energy of the coupled electron-lattice system E{ui...U]^; g) can be calculated from the 
equation: 

dE{ui...UN)/dg = ^Y^Uin{Ri) (3) 

i 

n{Ri) = {ip'^{Ri)iJj{Ri)) is the number of electrons per unit cell at the point R^. In terms of 
the electronic Green function (at T = 0, (jj), GapiR^t; R' ,t') = -i{T{^paiRi,t)^^{R',t')), 
one has: n{Ri) = —2iG{R,t; R,t + 6). (The Green function is diagonal in the spin indices: 
G api^R^t] R\t') = — i6ai3G(^R,t] R' ,t') ). 



III. COUPLED ELECTRON- LATTICE SYSTEM IN THE CLASSICAL REGIME. 

At static displacements, one needs to know only the frequency component, G{R, R';co): 

/I 
— G{R,R']u)exp{-iuj{t-t')) (4) 
27r 



Re-arrange the power expansion of G{R, R'; u) in Ui as: 



G{R, R!- u) = Go{R - R'; w) + ^ G'o(i? - Rf, uj)gu^Go{R - Ru u) (5) 

i 

+ ^ GoiR - Ru uj)guiGo{Ri - Rk, uj)gukGoiRk - R'; uj) + ■■■ 

The e-ph contributions in all powers in g are now summed first for the single site: Ui ^ Ui. 
Each of Ui,Uk--- in Eq.([5]) stands for the "dressed" local deformation at the corresponding 
site, Ri, Rki--i cind is determined by the relation: 

Ui=Ui + UiGo{Ri = Ri] u)gui (6) 



The free Green function, Go(-Ri = -Ri;i^)ll2| equals: 



Go{Ri = Ri] uj) = j {u - E{p) + Ef + tstgn{u)6)-' (7) 
With the help of identity 

(w — ^ + isign{uj)6)~^ = P{ } — iTxsign{uj)5{uj — ^) 

OJ — ^ 

and assuming the electron-hole symmetry for |^| = \E{p) — Ep\, one obtains: 

Go{Ri = Ri'ju) = —i7ih'{Ep)sign{ijj) (8) 

After trivial calculations, one finds the electronic contribution into the local elastic matrix 
at the site Ri, E{ui) = —2i dgUiG{Ri = Ri] t = t' + 5] g): 

E{u,) = -{W/27f) ln{l + {7ru{EF)gu,Y} (9) 

{W is the bandwidth). E{ui) is determined as the energy per one ion, v{Ep) is the number 
of states at the Fermi level per unit cell: ^{Ef) = ) ^o, where Sq is the area of the 2D 
unit cell. WulEp) = s is an insignificant model parameter, s ~ 1. We take s = 2. 
The interaction between two sites, E{ui, Uk), is given by the second term in Eq.Q: 

/d(jj 
— Gl{Ri^k]Oj) (10) 

{Ri,k = \Ri — Rk\)- 

The analytic form of G^lRi^k] ^) being cumbersome in 2D, for the estimate we use its 
asymptotic at PFRi,k > 1: 



/ 2tt 

Go{R]Uj) ^iu{EF)J^^ exp{isign{u)[{pF + {uj/vf))R - (11) 

to obtain: 

, 2 2/7-. \/ ,sin(2pp'i?) , , 

E{ui,Uk) = UiUug V {Ef){pfVf) (12) 

The total elastic energy of a single ion is the sum of Eqs. (P) and ([9]). Introducing: 

,^=f^^V^ (13) 



{MujIv{Ef)) = - (14) 



and the dimensionless notations for the ion's displacement, Ui = {ui/uo), one writes down 
the local elastic energy in the following simple form: 

Utotiui) = {l/7riy{EF)){u^, - ln[l + A^^]} (15) 

There, is the square of the dimensionless e-ph coupling constant and mq determines 
the spatial scale of the local elastic potential. The energy scale T* = I/ttu^Ef) = W/2ti 
is expressed through the electronic DOS. For 2if— NbSe2 the band calculations [l^ gave 
^{Ef) = 2.8 states /eV. Hence, T* is of the order of tenths of 1 eV. 



At A^ > 1 the potential Utot{ui) has two deep minima at u^- = ±y/l — A ^: 

Ut,tiu+,.) = {l/nu{EF)){l - A-2 - In A2} (16) 
At A^ < 1 and temperatures below T* -Uj ^ 1 and Utot{ui) can be written as: 

= {i/MEF))mi - A^) + Y (17) 

The quartic term in Eq. (|T711 is small, but the anharmonic contribution into the elastic 
ener gy i s the necessary ingredient in the molecular field approach to a CDW transition (see 



e.g. 



15|) 



Re-writing(pirfi?) as: (pfVf) = {^/'^^{Ef)){p^f^qI'^)i where 5*0 = (A/3/4)a^ is the area of 
the triangular unit cell, Eq. ( [T2l) in the same notations of Eqs. ( |T3|T4l) takes the form: 
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The numerical factor in this expression shows that the inter-site interactions are weak com- 
pared to the on-site Utot{ui), so that the structural transformation, if any, will occur at 
TcDW well below T*. At > 1 the model reduces to the model of interacting Ising spins: 



U- 



Vl - A-V„ (a, = ±1). 



When A^ < 1, the expressions for Utotiui) and Eq. ( 1T8|) for E{ui,Uk) together complete 
formulation of the problem: in the classical regime the system is fully described by the 
partition function Z(T,g): 

Z{T, 9) = j (ndui) exp[-^ ^ Utotiui, Uk)] (19) 

i,k 

( Utotal{Ui,Uj) = UtotiUi) + E{Ui,Uk)). 

The phase transitions are commonly treated in the molecular field approximation (see, 
e.g., [l5|). The method is not exact for short-ranged interactions, such as E{ui,Uk) in Eq. 
ffTSjl and we do not stay on details of the CDW transitions in TM2HDC. 

For the local properties one has for the partition function, Z{T,g): 

Z,{T,g) = j d^,exp[-— ^^-^K2-ln(l + A2s^)}] (20) 

Scattering of electrons on the lattice displacements above Tcdw is characterized by the 
average (gui) that enters the denominator of the Green function: 

G~\R; uj)=uj- E{p) + Ef + {gUi) (21) 

From Eq. (E]) it follows: gUi = gui{l — isign{uj)Ttv{Ep)gUi\^'^ . In the normal phase terms 
odd in Ui can be omitted. In dimensionless variables: 

/ 1 \ A^M^ 

gu. = ^s^gn{u) J (22) 

At A^ < 1 the imaginary part in G^^{R]u) is (gUi) = isign{uj)h?{'K /2)v{Ep)T: in the 
classical regime above Tcdw the resistivity of the system would be linear in T. In the 
opposite limit of A^ > 1, -u^ = -u^ _ = 1 — A^^, and the imaginary part is a constant of order 
of ~ l/'Kv{Ep) . The entropy for the system of non-interacting Ising spins is finite. 

In the ordered state with all ions occupying same minima, non-zero gUi 7^ stands 
together with the chemical potential by: 

gUi = ±(l/7rz/(E^))VA2 - 1 (23) 



7 



IV. DISCUSSION OF THE RESULTS 



A. Classical regime 

It is interesting to discuss relevance of the above to experimental results, in particular, 
for 2H-TaSe2. 

First, recall that the non-linear potential (Q was derived for a single ion. For Eq.(|9]) 
to be meaningful, uq must be small compared to the lattice parameter, a ~ 3.45A. With 
Wo ~ lOmeV taken for niobium, Eq. (fT^ gives Uq < 0.5 A. This justifies the assumption. 
Experimentally, shifts of cations in the superlattice lie in the interval between O.lA [t] and 
0.5 A [91, i.e. have same order of magnitude it should be for the model of Ising 

spins. 

The CDW-gaps seen in the tunneling non-linear current characteristics, dl /dV [g] (at T = 
A.2K) were large: ~ dOmeV and ~ 60meV for 2iy-TaSe2 and 2if-NbSe2, respectively. Such 
large gaps admit the interpretation in terms of energy of the bound electron-ion complex. 
As mentioned above, the energy scale, T* = 1/t:u{Ef) lies between 0.1-0.3 eV, depending on 
the number o: 
per unit cell 
value at > 1. 

A broad structure in the infrared conductivity with a threshold at ~ 0.25eV was discern- 
able up to 80K [5] (2i/-TaSe2). In recent ARPES experiments [161] features related to the 
CDW gap in 2if-NbSe2 were observable even in the normal phase at 115K. 

Above TcDW the imaginary part in Eqs. (1^ 1221) is large~ \/tiu{Ef), i.e., formation 
of the local electron-ion complex involves electrons from the considerable part of the band. 
At T < TcDW the CDW order parameter is proportional to ionic displacements in the 
superlattice oc gu{Q). No nesting is involved and, hence, no energy gap appears at the 
Fermi surface. 

(In the CDW-ground state constructed in [itI for the tight-binding model with the pre- 
vailing hopping matrix elements between the next- nearest neighbors, the FS at the F-point 



states per unit cell in the specific material (for NbSe2 v{Ep) ~ 2.8 States /eV 
l^).The value of \Utot{u+,-)\ may be adjusted even better to the observable 



17[ the CDW transition itself 



of BZ would remain metallic either. Note, however, that in 
takes place due to nesting on other parts of the electronic spectrum). 

The frequency of the collective excitations behaves differently at T = Tqdw for the local 
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potentials with one or two minima 16| (see also in 15|). In the former case the frequency 
vanishes at the temperature of the transition. In the latter, the frequency remains finite as 
an ion is now "trapped" by one of the two minimum. Neutron experiments jl|(b) for 2H- 
TaSe2 gave uj'^iQ) — 2Q{meVY at Tcdw- Softening of the acoustic phonons in 2i/-NbSe2 at 



TcDW = 33.57^ was observed, instead, in 



The difference in the phonon modes' behavior in the two materials needs a clarification. 
Recall that the longitudinal acoustic branch is coupled linearly with the optical mode, u 
of the same Si- symmetry at Q = a*/3. The Free Energy then has a contribution of the 
form:F(u, s) = u^ls"^ /2)+ul^{u'^ /2)+tus, where s stands for the acoustic branch. Minimizing 



F{u, s) and excluding u gives F{u, s) 



to: 



(t^/cj^)](s^/2). Let the optical mode, u be 



the mode that drives the transition. At ^^(T) (one minimum), the acoustic mode is 
the first one that manifests onset of the transition. If the potential has a few minima, the 



effective frequency of the acoustic mode w^^e// 



[w^ — (t'^/ul)] at Tcdw may or may not 
be zero depending on the temperature behavior of w^(r). Data [l](b) for 2H- TaSe2 agree 
better with the second possibility. 



The two 2if— material have different masses of Ta- and Nb- ions. It is known 16| that 
at lower temperatures, when quantum effects become important, cj^(T) — > even for a 
potential with a few minima. This would make possible the vanishing of the acoustic phonon 
frequency at Tcdw = 33. Si^' in the 2if-NbSe2 PJj and the finite oo'^iQ) at Tcdw = ^22K 



in 2iJ-TaSe2 



l|(b). 



B. Quantum regime 

At low T quantum effects become important in few aspects. Consider first the Schrodinger 
equation for an ion moving in the rigid potential Utotiui) = (l/7rz/(i?i?))f7(ttj), with U{ui) in 
the dimensionless notations. With the kinetic energy in the same notations, one has: 

~ i ^ ^(^) + (Uiu) - E)^{u) = (24) 

The dimensionless "mass", M in f l2^ is defined by 1/M = {'ki'{Ef)ujo)'^ /2. (The adiabatic 
parameter I/Vm = (tti'(Ep)ujq) / \/2 is about one tenth in 2/7— NbSe2). 

For the single-minimum well (A^ < 1) quantum effects are important when Tcdw < ^o- 
The lattice oscillations are now quantized and the electronic Green function is "dressed" by 
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phonons. In the adiabatic approximations, solution for the problem of interacting electrons 



18|- 



and phonons in normal metals was given many years ago in 

Eq. f l24p adds a possibility of quantum tunneling between minima of the two-well poten- 
tial. However, Eq.f l24p can only describe tunneling in a "rigid" potential, a potential built 
only by the interatomic forces in the lattice. Meanwhile, the two-well potential Utot{ui) in 
Eqs.l pfTSj) comes about selfconsistently due to interactions between electrons and the local 
lattice distortions. In a tunneling event, at which the ion goes, say, from m+ to the 
electronic configuration reverses as well. Such a feature cannot be described in terms of 
Eq.f l2^ . Besides, in the quantum regime Utgt^{ui) itself varies with temperature js]. 



The collective tunneling was discussed in [9| in connection with the martensitic transitions 
in NbsSn and VaSi. The authors were concerned with its mapping on the Kondo problem. 
(Another analogy would be the two- level quantum system). No closed solution was obtained 
in {9 1 for an effective tunneling matrix element. 

Below the transition in the CDW ground state ions will occupy their proper minima, in 
accordance with the superlattice pattern [Ll](b) and tunneling between minima must stop. 



Nevertheless, experimentally, the pattern reveals an interesting peculiarity. Indeed, in the 
CDW phase one in three atoms along the symmetry lines finds itself in the position with the 
trigonal symmetry (2H- TaSe2 |l](b); see Fig.4 in js]). While below transition the intersite 
interactions do indeed arrest quantum tunneling between minima for the other two of the 
three atoms, the degeneracy is not lifted for the atom in this symmetric position and the 
collective quantum tunneling processes between minima continue at this site down to T = 0. 

According to ARPES data [19], the CDW transition in 2H- NbSe2 slightly affects only 
the two-barrel Fermi surface at the K-point in BZ, while FS at the F - point remains intact. 
(Three Q-vectors couple together the three points on the inner FS at K, which is seen 
experimentally as small but observable gaps at these three points. The small local gaps 
(~ 2AmeV) survive even in the normal phase at 115K. This result is not quite clear.) 

The physics of strong local e-ph interaction bring us back to the Holstein model [l3| of 
electrons interacting with dispersionless phonons. No exact solution is known for the Hol- 
stein model either, but its low teraperature physics was investigated numerically in DMFT 
(Dynamical Mean Field Theory J20|) approximation. (For a brief summary of results for 
the Holstein model at T=0 see |2l|). 

The DMFT approach is strictly local. Intersite interactions (see Eq.f llSp ) and the CDW 
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transition itself cannot be treated by DMFT. 

Here we indicate parallels between our physics above and the results [21I] at T=0 for the 
Holstein model. Among them: 

1) The double- well potential develops when the e-ph coupling constant becomes larger 
some critical value; 2) Large imaginary part in the Green function: the Spectral Function 
extends over the energy interval that significantly exceeds the phonon frequency (compare 
with our Eq. f l22|) ): 3) The ground state remains metallic (at least for not-too-strong e-ph 
coupling^ . 

In 2l| the mass of electronic excitations increases as the residue at the pole of the 
electronic Green function decreases. Judging by these results, the physics studied above 
allows its extension into the low temperature regime. 



V. SUMMARY 



In summary, strong e-ph interactions may change the local environment in the lattice 
qualitatively by binding ions and electrons together, thus breaking the adiabatic approx- 
imation [18]. In such case the CDW transition is the phase transition in the system of 
interacting Ising spins. The transition is the result of intersite interactions that come about 
due to the exchange by electrons between local sites. The sizes of the structural vector Q 
and FS being not related to each other, the energy spectrum of electrons remains metallic 
below the CDW transition. Large values of the CDW gaps observed experimentally receive 
the natural explanation in terms of deep energy minima of the two-well local potential. 
Conclusions from the analytical results derived in the classical regime can be extended to 
lower temperatures as it follows from the comparison with numerical results for the Holstein 
model. 
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